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ABSTRACT 

HD 177830 is an evolved K0IV star with two known exoplanets. In addition to the planetary 
companions it has a late-type stellar companion discovered with adaptive optics imagery. We observed 
the binary star system with the PHARO near-IR camera and the Project 1640 coronagraph. Using 
the Project 1640 coronagraph and integral field spectrograph we extracted a spectrum of the stellar 
companion. This allowed us to determine that the spectral type of the stellar companion is a M4±1V. 

We used both instruments to measure the astrometry of the binary system. Combining these data 
with published data, we determined that the binary star has a likely period of approximately 800 
years with a semi-major axis of 100-200 AU. This implies that the stellar companion has had little or 
no impact on the dynamics of the exoplanets. The astrometry of the system should continue to be 
monitored, but due to the slow nature of the system, observations can be made once every 5-10 years. 

Subject headings: binaries: visual - instrumentation: adaptive optics - stars: individual(HD 177830) 


1. INTRODUCTION 

HD 177830 (HIP 93 746 = W PS 19053+2555) is an 
evolved K0IV star ( Vogt et ahl 120001 1 that is on the 
boun dary between giants and subgiants (|Ghezzi et al.1 
HM3). It i s a nearby star with a distance of 
59.0+2.2s pc l|van Leeuwenl I2007H . There have been 
several age determinations including an age of 3.8l^g 
Gyr (j Valenti fe Fischerl 12005H . 4.03 Gyr dSaffe et al.l 
[2(T05h 3,28 f A Gvr dTakeda et al.ll2007lh and 3.46+0.29 
Gvr (Uo fre et, a ll 1201511 all of which are consistent. 
iBaines et al.l (l2008al ) measured the diameter of the star 
interferometrically with the CHARA Array as 2.99+0.39 
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R© and estimated th e T e f f as 4804 K. Using high resolu- 
tions visible spectra, iMortier et al.1 (I2013H measured the 
metallicity, log g 7 T e ff, and then used these param¬ 
eters to estimate the mass using t heoretical iso chrones 
and a Bayesian estimation method (Ide Silva et al.l[2006ii . 
This resulted i n a m ass of 1.17+0.10 Mq, and T e ff 
as 4752+79 K. Uofre et all (|2015lf carried out a sim¬ 
ilar analysis using different input spectra and model¬ 
ing tools to derive the input parameters, but used the 
same software to compute masses. Their estimates of 
1.37+0.04 M 0 and T e// of 5058+35 K are slightly dif¬ 
ferent. The system has been observed with multiple in¬ 
struments on the Spitzer telescope and none of the ob- 
servations have turned up any evidence o f a debris disk 
(iBeichman et alJl2005l iTrilling et aLll2008t fBrvden et al.1 
l2009UTanner et alJl2009t iDodson -Robins on et al.ll2011h . 

By mo nitoring radial velocity (RV) variations, 
IVogt et all (|2000f ) discovered a planet orbiting 
HD 177830 with a period of 391.6 days, an eccen¬ 
tricity of 0.41 and an M sin i of 1.22M,/, where M is 
the mass and i is the inclination. This planet was 
designated HD 177830b. It orbits the primary at a 
distance of 0.63-1.57 AU, and lies inside the h abitab le 
zone. The planetary or bit was updat e d bv | Butler et akl 
(|2006H and again by iWright et al.l (|2007h who also 
showed the first signs of a second plane t in the system. 
The s tar continued to be observed and iMcschiar i et al.l 
(j2011f) published an update of the orbital parameters 
of HD 177830b and announced the discovery of an 
inner planet in the system. This planet, HD 177830c, 
has a period of 110.9+0.1 days and eccentricity of 
0.3+0.1 and an Msini of 0.15+0.02 Mj. The updated 
parameters of HD 177830b are a period of 407.31 days, 
an eccentricity of 0 .009+0.004 and an M sini of 1.48Mj. 
[Reffert fc Quirrenbacbl (j2011l) were able to use the lack 
of astrometric detection with Hipparcos satellite to place 
an upper limit on the mass of HD 177830b of 225.2 Mj. 
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lEggenbereer et all (12007ft first reported the discovery 
of a stellar companion to HD 177830 with a separation of 
1"6 using adaptive optics (AO) on the Very Large Tele¬ 
scope. The stellar companion was designated HD 177830 
B, while the primary was designated HD 177830 A. 
Based on near-IR colors and the differential magnitude 
of 6.6 at 1.6 (im they estimated the spectral type of 
HD 177830 B to be M2V-M 5V with a mass of 0.23±0.01 
Mg. iRoberts et al.1 (l201lft published an observation of 
th e same companion t hat pr edates the discovery images 
of lEggenberger et al.1 (|2007jl using the visible-light AO 
system on the Advanced Electro-Optical System (AEOS) 
telescope. The astrometry from both papers is shown in 
Tabic m 

There have been several prior high angular resolution 
observations of the HD 177830. iBaines et al.1 (I2008bft 
used the CHARA Array to search for stellar companions 
and did not detect the stellar companion because it was 
outside of the instrument’s field of view of O'/I. Inside 
of their field of view, they were able to rule ou t any ad- 
ditional stellar companions earlier than K0V. iLu et all 
(11987ft carried out speckle interferometry observations 
with the Kitt Peak 4 m telescope, but the observations 
did not achieve a sufficient dynamic range to detect the 
companion. 

The stellar companion’s measured separation corre¬ 
sponds to a projected separation of 97AU and this sep¬ 
aration has the potential to have an impact on the 
dynamical behaviour and evolution of the exoplanets 
orbiting HD 177830 A. The population of exoplan¬ 
ets in binaries with semimajor axes smaller than 100 
AU is st atistically different than those orbiting sin¬ 
gle s tars dZucker fe Mazehl l2002t iBonavita fe Desideral 
HMD, while for binaries with separations larger than 
100 AU the frequencies of exoplanets among single 
stars and components of wide binaries are indistinguish - 
able (|Raghavan et aljr2006l : IBonavita fe Desideral 12007ft . 
With that in mind, we carried out observations with two 
instruments on the Palomar 5m telescope in an effort 
to understand the orbital parameters of this exoplanet 
hosting binary star. 

2. OBSERVATIONS AND DATA ANALYSIS 
2 . 1 . Project 16j0 

HD 177830 was o bserved with the Project 1640 
(P1640) coronagraph (jHinklev et al.lf20TTI l on 2012 June 
12 UT. The P1640 coronagraph is integrated with an in¬ 
tegral field spectrograph covering the YJH bands. Like 
PHARO, this instrument is mounted on the PALM-3000 
AO system. We collected 3 data cubes, each with an inte¬ 
gration time of 549.9 s. For each data cube the primary 
was placed behind the occulting disk. The data were 
reduced using the Projec t 1640 data reduction pipeline 
(Zimmerman ct al.1 12011ft . Since that paper was writ¬ 
ten, a few upgrades have been made to the pipeline that 
reduce lenslet-lenslet cross talk. The pipeline processing 
produces an image of the object in each of 32 wave bands 
resulting in a data cube. 

While we collected several unocculted images of 
HD 177830, the dynamic range between the companion 
and primary makes it extremely difficult to detect the 
companion. Two of our occ u lted ima ges used the as- 
trometric grid spots dSivaramakrishnan &; Onnenheimed 


12006ft . We fit a Gaussian to each of these spots. Then 
we identified the location of the primary by fitting lines 
to the vertical and the horizontal spots. The primary is 
located where these two lines intersect. Figure Q] shows 
an image slice with the grid spots. The final astrometry 
is the average of the astrometry measured from 26 of the 
data slices from each of the two data cubes with grid 
spots. The frames in the water absorbtion bands with 
center wavelengths of 1145 ran, 1170 nm, 1395 nm, 1420 
nm, 1445 nm and 1795 nm had too low signal to fit a 
Gaussian to the grid spots. The astrometry is presented 
in Table [lj which lists the Besselian date of the obser¬ 
vation, the position angle ( 9 ) and separation (p) of the 
system and the instrument used in the observation. 



Figure 1. A slice of the occulted P1640 image cube at a wave¬ 
length of 1.495/rm. It has been rotated so that North is up and 
East is to the right. There are four astrometric grid spots in the 
image, see the text for a discussion of their use. An arrow points 
towards the companion. The field of view is 4 '.' 0 . 

2.2. PHARO 

We observed HD 177830 on 2012 May 9 UT and on 
2014 May 14 UT with the Palomar Observatory Hale 5 
m telescope using the PALM-3000 AO system and the 
PHARO near-IR camera. The PALM-3000 AO system 
is a natural gu ide star sy stem using two deformable mir¬ 
rors (DM) (|Dekanv et alJlMToli . One DM corrects low- 
amplitude high spatial frequency aberrations, while the 
other corrects the higher-amplitude low spatial frequency 
aberrations. The system is optimized for high contrast 
observations and routinely produces Strehl ratios greater 
than 80% in the Ks band. The PHARO camera uses 
a HgCdTe HAWAII det ector for observation s between 1 
and 2.5 pm wavelength (Havward et al.lfcOOlft . The cam¬ 
era has multiple filters in two filter wheels and we used 
the Ks filter to collect 10 frames in 2012 and 50 frames 
in 2014. In 2012, only sky frames were collected for cal¬ 
ibration purposes. Past experience with PHARO has 
shown that its flat fields and dark frames are very stable 
from run to run. We used calibration data from 2013 to 
calibrate the 2012 images. Calibration frames were col¬ 
lected on the same night as the 2014 data. The night of 
2014 May 14 UT suffered from high winds resulting in 
poor seeing and lower image quality than is normal for 
P3K. After debiasing, flat fielding, bad pixel correction 
and background subtraction, the frames were coadded. 
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The resulting images are shown in Figure [2] The fit- 
stars algorithm was used to measure t he astro metry and 
photometry of the objects (Iten Brummelaar et al.l fl996l 
1200011 . Photometric er ror bars were assign ed using the 
technique described in lRoberts et al.1 (12005H . Astromet¬ 
ric error bars were set equal to the median error bar of 26 
binaries measured with P3K and PHARO ([Roberts et al.1 
I2015H . The resulting measurements are in Table |T) It has 
the same format as the P1640 results, but we also list the 
differential magnitude between the two stars in the Ks 
filter. 



Figure 2. Ks images of HD 177830 from the PHARO instrument 
taken in 2014. An arrow points to the companion. In the image, 
North is up and East is to the right. The close object to the south 
east of the primary is a ghost from a neutral density filter in the 
camera. This is an approximately 4" wide subimage from the full 
25"field of view. 


3. SPECTRAL TYPE DETERMINATION OF THE 
COMPANION 

To create a spectrum of the stellar companion, we per¬ 
formed aperture photometry on the images with suf¬ 
ficient S/N in the P1640 data cube. We used the 
same tech niques that we have used in our previous pa¬ 
pers (e.g. iZimmerman et al.1 120101 : iHinklev et al.l l2010t 
IRoberts et all 120121 ). This was done with the aper.pro 
routine which is part of the IDL astr olilP^l and is an 
adaptation of DAOphot (iStetsonl 11987T) . A photometry 
aperture and a sky annulus were centered on the com¬ 
panion. The radius of the photometry aperture was set 
equal to 3.5A /D rounded to the nearest pixel size, where 
A is the central wavelength of each image and D is the 
aperture of the telescope. This corresponds to a radius 
of O'/18 at the center of J band and 0'/24 at the center 
of H band. The radius of the sky anullus was set to 
0.67"; this is large enough to avoid significant portions 
of the central PSF while also avoiding the occultation 
spot. The width of the sky annulus was set to 10 pixels, 
or 0.19". The background was set equal to the average 
intensity of the pixels in the sky annulus and was sub¬ 
tracted from each pixel in the photometry aperture. The 
spectrum is the summed power in each slice as function of 
wavelength. The aper.pro code also produces a measure¬ 
ment error bar based on photon statistics and the error 

14 http://idlastro.gsfc.nasa.gov 


in the background estimation. The extracted spectrum 
is a convolution of the object spectrum and the spectral 
response function (SRF). The SRF includes the instru¬ 
mental response and the absorption due to the Earth’s 
atmosphere. We computed the SRF by using an unoc¬ 
culted observation of HD 177830 A as a reference source. 
We extracted the spectrum of HD 177830 A using the 
aperture photometry method described above. Then we 
divided the extracted spectrum by a reference spectrum 
of the same spectral type. HD 177830 A is a K0IV. We 
did not have a template for this spectral type, instead 
we used a K0V spectra from the In frared Telesco pe Fa¬ 
cility (IRTF) Spectral Library (|Ravner et al.ll2009l ). This 
produces an estimate of the SRF. The companion’s mea¬ 
sured spectrum is then divided by the SRF and yields 
the calibrated spectrum of the companion. The portions 
of the spectrum that overlap with the water bands were 
removed, as they had large errors. 

Since the template spectrum had the incorrect lumi¬ 
nosity class, we checked to see how much impact the 
incorrect luminosity class had on the resulting spectra. 
We computed the spectrum of HD 177830 B with an SRF 
from HD 129814 (spectral type G5V), taken two nights 
later at a similar air mass. The resulting spectrum for 
HD 177830 B is the same to within the error bars. 

The error bars on the final reduced spectra are the 
combination of the errors of the measured science spectra 
and the error on the SRF. The SRF errors come from the 
combination of the errors of the measured calibration 
spectra and the errors on the template spectrum. See 
IRoberts et al.1 (1201211 for further discussion on SRFs. 

After the spectrum of the companion was extracted 
from the P1640 data, we comp ared it aga inst the spectra 
in the IRTF Spectral Lib rary dCushing. R.avner fc Vaccal 
120051 : iRavner et all 20091 1. These include FGKM main se¬ 
quence stars and LT brown dwarfs. The template spec¬ 
tra were binned and smoothed in order to produce the 
equivalent spectra to having the star observed by P1640. 
Each spectrum was normalized by its average. Then each 
template spectrum was compared against the measured 
spectrum using the sum of the squares of the residual 
(SSR) as a metric, 

SSR = J2 ™\(S x ~Rx) 2 , ( 1 ) 

A 

where, S\ is the measured spectrum at a given wave¬ 
length, R\ is the binned reference spectrum at the same 
wavelength and w\ is the weight at the wavelength. The 
best fit reference spectrum was the one with the mini¬ 
mum value to the metric. The weights were set equal 
to the inverse of the computed error at each wavelength 
point. 

The extracted spectrum of HD 177830 B is shown in 
Figure [3] The error bars for each data point are shown 
in the figure. The error bars are smallest in the 77-band 
portion of the spectrum and since the weight of each data 
point is the inverse of the error bar, this means that the 
77-band data are weighted more heavily than the J-band 
data. This is appropriate for several reasons. The AO 
performance is improved in the 77-band and produces a 
higher quality image. The detector samples the point 
spread function (PSF) at a higher spatial frequency in 
77-band and finally the coronagraph is optimized for 77- 
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Table 1 

Measured astrometry and photometry for HD 177830. 


Epoch 

e (°) 

pH 

A M ( Ks) 

Instrument/Reference 

2002.5474 

84.1±1.0 

1.62i0.01 


Roberts et al. (2011) 

2004.4784 

84.85±0.21 

1.645i0.008 


Eggenberger et al. (2007) 

2005.3494 

84.60i0.39 

1.64i0.01 


Eggenberger et al. (2007) 

2012.3549 

84.3i0.2 

1.67i0.01 

6.5i0.3 

PHARO 

2012.4475 

86.0i0.1 

1.68i0.002 


PI 640 

2014.6748 

85.3i0.2 

1.67i0.01 

6.1i0.3 

PHARO 

band observations. This results in 

a slightly higher con- 

ondary corresponds to a main sequence mass of 0.2 M fl 

trast in H-band (IHinklev et alJ 2011) and lower errors. 

(Reid & Hawlev 2005|). 

In Figure [3] we have over plotted the spectra for M2V, 

As discussed in the introduction, there are two lit- 


M4V, and M6V stars. The best fit (via Equation [T]) for 
the companion’s spectral type is M4V. FigureS] shows the 
SSR metric plotted as a function of spectral type. Several 
adjacent spectral types have similar fitting metric values 
and their spectra in Figure[3]appear to fit the data almost 
as well within the error bars. From this we conclude that 
the spectral type is M4±1V with a conservative error bar. 
This agrees with t he earlier photometric det ermination 
of M2V-M5V from lEggenberger et al.l (120071) . 

4. ORBITAL ANALYSIS 

lEggenberger et all ( 12007? ) concluded that the binary 
star system had co mmo n prop er motion. With the ad¬ 
ditional data from Roberts et al.l (120121 ) and this pa¬ 
per, that conclusion is strengthened. Between 2002 and 
2014 the companion moved at a rate of 4.8 mas/yr rel¬ 
ative to the primary, while if it was a fixed background 
star it w ould have been expected to move -65.92 mas/yr 
ijvan Leeuwenll2007f) due to the proper motion of the pri¬ 
mary. This clearly shows that the companion has com¬ 
mon proper motion with the primary. We plot the rela¬ 
tive motion in Figure[5] From this we can tell the bulk of 
the motion is to the East, but due to the size of the error 
bars in relationship with the small amount of motion, 
we can not tell if the companion has an increasing or 
decreasing declination. The 2012.4475 data point from 
P1640 does not line up with the rest of the data. This 
is probably due to an error in the plate scale for that in¬ 
strument. The separation corresponds to approximately 
90 pixels and a 1% error in the plate scale would throw 
the measurement off. 

Using the technique of iTokovininl (12014D . we can con¬ 
strain the period of the binary. Using the equation, 


P* = 


P 3 tt - 3 

Mi + M 2 


1/2 


( 2 ) 


erature estimates of the mass esti mates from spectro- 
scopic modeling: 1.17±0.10 M n (IMortier et ahl 120131) 
and 1.37± 0.04 (Uofre et al1l2015l) . Since the T e ff 
derived by IMortier et al.l (120131) matches the interfero- 
metrically derived T» o oflBaines et all (l2008al) . we used 
the mass estimate of IMortier et alJ (|2013li resulting in 
a peri od es timat e of 829 yrs. The mass estimate of 
Uofre et all (1201511 produces a period of 771 years. 

With a semi-major axis of 100-200 AU the stellar com¬ 
panion will not have influenced the formation or dynam¬ 
ical evolution of the exoplanets unless the stellar orbit 
is extremely eccentric. There is little evidence that the 
stellar companion has impacted the dynamics of the exo¬ 
planets, with semi-major axes of 1.2218±0.0008 AU and 
0.5137±0.0003 AU. Most of the suggested ways that bi¬ 
nary companions interact with exoplanets result in plan¬ 
ets with eccentric orbits dKlev fe Nelsonl I2008D or hot 
Jupiters (|Wu fe Murravl 120031 : INaoz et al.ll2012D . In this 
case, HD 177830 b has a low eccentricity of 0.009±0.004 
and appears not to have had its orbit altered by the stel¬ 
lar companion. The inner planet, HD 177830c, has a 
low to moderate eccentricity of 0.3, but since it is inside 
of the orbit of HD 177830 b, it seems unlikely that the 
stellar companion would have effected its orbit and not 
that of HD 177830 b. It is possible that the system has 
additional undetected planets in outer orbits that have 
been modified by the gravitational pull of HD 177830 
B. It is also possible that the system used to have ad¬ 
ditional outer planets that were ejecte d from the system 
due to interactions with the binary dKaib et al.l 1201 ill ). 
Since the eccentricity of HD 177830b is almost zero, the 
system does not appear to be undergoing a Kozai cycle, 
implying that the relative inclination of the stellar com¬ 
panion’s orbit and t hat of the p la nets is less than the 
critical angle of 39.2° (|Holman et~aLlll997D . 

5. NULL RESULTS 


where, P* is the probable period based on Kepler’s third 
law, p is the measured separation, assumed to be the 
semi-major axis, 7r is t he para llax, and M\ and M 2 are 
the mass of each star. ITokovininl (120141) carried out sim¬ 
ulations with random orbital phases and eccentricities 
and showed that the median ratio of p/a is close to 1 
depending slightly on the eccentricity distribution used 
in the simulations. In all cases, the ratio never exceeds 
2. This results in the ratio of P*/Ptrue always being less 
than 3.17. For the case of HD 177830, we use the latest 
separation measureme nt from Table [1 ] an d the Hippar- 
cos measured parallax (Ivan Leeuwenl 120071) . For masses, 
the spectral type of M4V that we derived for the sec- 


In addition to HD 177830, we used P1640 to observe 
several exoplanet host stars that did not have known 
companions. These observations happened during the 
first phase of the P1640 instrument w hen it was use d with 
PALM-3000 pre decessor, PAL AO (|Dekanv et alJ 119971 : 
iTrov et aHl2000if . We did not detect any candidate com¬ 
panions in our 3'/4 held of view. Table [2] lists the name 
most commonly used in reference to the exoplanet, the 
Hipparcos number, the Besselian date of the observations 
and the radius of the held of view in terms of AU. Af¬ 
ter data reduction with the versio n of the LOCI sp eckle 
suppression software described in dCrepp et al.ll201ll ). all 
the stars had contrast curves approximately the same as 
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Figure 3. The spectrum extracted from the P1640 data of HD 177830 B. The three over plotted template spectra are Green=M2V, 
Orange=M4V, and Red=M6V. The determination of the spectral type is M4±1V. 
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Figure 4. The SSR metric from Equation [T] plotted as a function 
of spectral type for HD 177830. The units of SSR are arbitrary. 
For some spectral types, the IRTF Spectral Library has multiple 
stars of that class; this results in multiple data points for several 
of spectral types. The best fit spectral type is M4V. 


Table 2 

Exoplanet Host Stars with No Detected Companions 


Name 

HIP 

Date (UT) 

FOV Radius (AU) 

HD 5319 

4287 

2008.8058 

195 

HD 69830 

40693 

2009.2054 

21.2 

GJ 849 

109358 

2008.5242 

15.5 


that shown in that paper. The contrast is high enough 
to rule out additional stellar companions within the field 
of view. 

6. SUMMARY 


Figure 5. A plot of the relative motion of the companion to the 
primary star. The primary direction of relative motion is to the 
East, but we can not tell if the declination is changing or not. 


We used the P1640 coronagraph and IFU to determine 
that the secondary in the exoplanet hosting binary sys¬ 
tem, HD 177830, is a M4±1V star. We also extracted the 
astrometry from this measurement as well as two oth¬ 
ers taken with the PHAR.O instrument. Combined with 
published astrometry, we are able to constrain the orbit 
to a period of approximately 800 years and a semi-major 
axis of 100-200 AU. This strongly suggests that the bi¬ 
nary system has had little to no impact on the known 
exoplanets in the system. The computed orbit of the 
system can be improved with additional astrometry, but 
due to the slow motion of the system, the observations 
can probably be made every 5-10 years until significant 
orbital motion is detected. 
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